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Abstract. 

In this article we construct zcta functions of quantum graphs using a contour 
integral technique based on the argument principle. We start by considering the special 
case of the star graph with Neumann matching conditions at the center of the star. We 
then extend the technique to allow any matching conditions at the center for which the 
Laplace operator is self-adjoint and finally obtain an expression for the zeta function of 
any graph with general vertex matching conditions. In the process it is convenient to 
work with new forms for the secular equation of a quantum graph that extend the well 
known secular equation of the Neumann star graph. In the second half of the article 
we apply the zeta function to obtain new results for the spectral determinant, vacuum 
energy and heat kernel coefficients of quantum graphs. These have all been topics of 
current research in their own right and in each case this unified approach significantly 
expands results in the literature. 
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1. Introduction 

Zeta functions associated with combinatorial graphs have been well studied. The 
Ihara zeta function [21] was interpreted in terms of a finite graph by Sunada |44j . 
for subsequent generahzations see e.g. |2ll [2l |43]. Such combinatorial graph zeta 
functions are defined by an Euler product over sets of primitive cycles (periodic orbits) 
which are itineraries of graph edges. A quantum graph identifies bonds (edges) of 
a combinatorial graph with closed intervals generating a metric graph. In addition a 
quantum graph has an operator acting on functions defined on the collection of intervals. 
This is typically a self-adjoint Hamiltonian operator acting on functions in a Hilbert 
space defined on the metric graph; in the context of inverse spectral problems see, e.g., 
[9]. Quantum graphs were introduced to model electrons in organic molecules, some 
more recent applications include areas of mesoscopic physics, quantum chaos, photonic 
crystals, superconductivity, the quantum Hall effect, microelectronics and the theory 
of waveguides, see |37l |39] for recent surveys of the applications of graph models. In 
such applications the spectrum of the quantum graph determines physical properties of 
the model. In this sense the quantum graph is not far removed from its combinatorial 
cousin where important graph theoretic properties like the connectivity of the graph are 
approached through the spectrum of a discrete Laplace operator on the graph. 

Spectral questions asked of quantum graphs are traditionally approached through 
a trace formula. Trace formulae express spectral functions like the density of states or 
heat kernel as sums over periodic orbits on the graph. They naturally acquire a similar 
flavor to the Ihara zeta function. The first graph trace formula was derived by Roth 
|40j . Kottos and Smilansky [35] 136] introduced a contour integral approach to the trace 
formula starting with a secular equation based on the scattering matrix of plane-waves 
on the graph. Solutions of the secular equation correspond to points in the spectrum of 
the quantum graph. (For the current state of the art of graph trace formulae see [8].) 

If the spectrum of the quantum graph is ^ Aq ^ Ai ^ . . . then formally the 
associated spectral zeta function C(s) is defined by 



where the prime indicates that zero modes, if present, are omitted from the summation. 
Using the trace formula such a zeta function can immediately be written as a sum over 
periodic orbits. However, adopting this approach the zeta function and the trace formula 
are closely related and the zeta function does not provide new insight. Instead we take 
a step back from the trace formula and derive the zeta function using a contour integral 
based on new secular equations for the graph. Applying a contour transformation we 
obtain an integral formulation for the zeta function based on the matching conditions 
at the graphs vertices directly in the general case. The program can be thought of as a 
generalization of the special case of a star graph with Neumann matching conditions at 
the vertices (defined in Section [3]). In this case there is a well known secular equation 



oo 




(1) 



i=0 
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whose solutions define the spectrum of the star [51 E], 

B 

^tanA;Lfe = 0. (2) 

b=i 

The sum is over the B arms of the star and is the length of the 6'th arm. If kj is a 
solution of (|2]) then \j = k"^ is an eigenvalue of the Laplace operator on the star. In such 
cases the simple structure of the secular equation produces more exphcit formulations 
of the zeta function. 

The zeta function of a particular infinite quantum graph whose spectrum is 
related to Diri chiefs divisor problem is studied in [15] and the zeta function of 
the Berry-Keating operator appears in [T3], however there was previously no general 
formulation of zeta functions of the Laplace operator on graphs. In addition many 
spectral properties of quantum graphs, that are subjects of active investigation 
PI3[HllI2[13l[ni[l9l|25l[33l|40l|46],canbe obtained in a straightforward way from the 
zeta function. In the second half of the article we apply our representations of the zeta 
function to derive the spectral determinant, vacuum energy and heat kernel asymptotics 
of quantum graphs. In all these cases the zeta function approach adopted here yields 
new results. In detail we obtain the following theorems. 

Theorem 1. For the Laplace operator on a graph whose vertex matching conditions are 
defined by a pair of matrices A andM, with AM^ = MA'^ and rank(A, B) = B, the spectral 
determinant is 

2^/(0) 



det'(-A) 



b 



where 



/(t) = det A - tB 



coth(tL) — csch(tL) 
— csch(tL) coth(tL) 



Civ is the coefficient of the leading order term t ^ in the asymptotic expansion of f{t) 
as t — 7- oo. 

Theorem 2. For the Laplace operator on a graph, under the conditions of TheoremUl 
the Casimir force on the bond (3 is 



24L2 Trio dLp 

provided the graph is generic: the poles of f{z) = f{—iz) are the whole of the set 
{rmr/Lh\m G Z, 6 = 1, . . . , B}. 

In this article we follow the scheme introduced by Kostrykin and Schrader [34j to 
classify matching conditions of self-adjoint Laplace operators on the graph, so matching 
conditions at the vertices are defined by a matrix equation Ai/j -|- Bi/>' = where t/? and 
i/j' are vectors of the values of a function ip and its outgoing derivatives at the ends of 
the bonds, see Section [2j Theorems [1] and [2] provide new general formulations of the 
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spectral determinant and vacuum energy for the Laplace operator on a graph with any 
vertex matching conditions illustrating the power of this unified approach. 

The article is structured as follows. In Section [2] we define the quantum graph 
models that we study. Section |3] introduces the contour integral approach used to 
evaluate the graph zeta functions with the simplest example of a star graph with 
Neumann matching conditions at the center. Section H] generalizes the star graph results 
by allowing any form of matching conditions at the center consistent with a self-adjoint 
realization of the Laplace operator on the star. Having demonstrated the techniques used 
to incorporate general matching conditions at a vertex we formulate the zeta function of 
general quantum graphs in terms of the vertex matching conditions in Section |5] which 
is the central result of the paper. As a byproduct we also formulate a number of new 
forms for the secular equation of a quantum graph. 

In the second half of the paper we investigate the implications of our results 
for a number of spectral quantities. Section [6] compares the results for the spectral 
determinant derived from the zeta function to those already obtained for quantum 
graphs. The zeta function approach developed here is not only more direct but the 
results obtained have a particularly simple form being expressed directly in terms of 
matching conditions. Section [7] uses the zeta functions to obtain a new formulation of 
graph vacuum energy. Section [8] determines implications of the zeta function results for 
the asymptotics of the heat kernel. In Section [9] we provide a concrete and current 
application of our results to a graph theoretic generalization of a piston. In the 
Conclusions we point out the most important results of our contribution. 

2. Quantum graph model 

For the purpose of this article the particular quantum graph model we consider is a 
self-adjoint Laplace operator on a metric graph. However, the term quantum graph 
is often applied more widely to describe self-adjoint differential or pseudo-differential 
operators on metric graphs. In this section we introduce the graph models we employ, 
for a general review of analysis on quantum graphs see [201 [38] . 

A graph is a set of vertices connected by bonds, see for example Figure [TJ In a 
metric graph G each bond h is associated with an interval [0, L^] so Lh is the length 
of h. For a bond h = {v, w) connecting vertices v and w the choice of orientation for 
the coordinate Xb on the interval [0, L^] is arbitrary, our results are independent of this 
choice of orientation. However, for the sake of clarity, when a bond is written as a pair 
of vertices b = {v,w) the coordinate x;, = at t> and Xb = Lb at w. The vertices and 
bonds are enumerated so G {1, 2, . . . , V} and b G {1, 2, . . . , B}. The total length of 
the graph G is denoted by £ = '^b=i 

A quantum graph consists of a metric graph with a self-adjoint differential operator 
on the set of intervals associated with the graph bonds. In this article we consider 
Laplace operators on metric graphs. The differential operator on the bonds of the 
graph is — A function on G is defined by the set of functions {'ipb{xb)}b=i,...,B on 
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the intervals associated to the bonds. The Hilbert space for the graph is consequently 

B 

H = 0L2([O,L,]). (3) 

6=1 

A self-adjoint realization of the Laplace operator on G is determined by specifying a 
suitable domain in "H. This can be achieved by defining an appropriate set of matching 
conditions at the vertices of G. We assume all matching conditions are local: at a 
vertex v matching conditions respect the connectivity of the graph only relating values 
of the function and its derivatives at the ends of the intervals connected at v. General 
matching conditions at all vertices of G are specified by a pair of 2B x 2B matrices A 
and B. For a function %p on G let 

=(V'i(0),...,^B(0),V'i(i^i),...,^B(i^B))^, (4) 

V'' = (V^'i(O), . . . , M^), -V^ULi), . . . , -MLb))^ ■ (5) 
The graphs matching conditions are then defined by the matrix equation 

AV' + BV' = . (6) 

The following theorem of Kostrykin and Schrader classifies all matching conditions 
of self-adjoint realizations of the Laplace operator [3l]. (An alternative unique 
classification scheme was introduced by Kuchment in |38j.) 

Theorem 3. The Laplace operator with matching conditions specified by A and B is 
self-adjoint if and only if (A, B) has maximal rank and AM^ = MA^ . 

Although we set out a general classification scheme for self-adjoint Laplace operators 
on graphs it will not be used until Section HI initially we concentrate on the simpler case 
of a star graph with Neumann like matching conditions at the center of the star. 

Given a self-adjoint Laplace operator — A on a graph we are interested in properties 
of the spectrum Aq ^ Ai ^ A2 ^ . . .. It is convenient to introduce an alternative spectral 
parameter k so we study solutions of the eigenproblem 

- AV^ = k'^tp . (7) 

Then Xj = A;| and we refer to the non-negative sequence ^ /cq ^ /^i ^ A;2 ^ • • • as the 
/c-spectrum. Formally the spectral zeta function of a quantum graph that will be our 
object of study is, 

00 

j=0 

3. Zeta functions of star graphs with Neumann matching conditions at the 
center 



To demonstrate the technique we start with a model problem of a star graph with 
Neumann like matching conditions at the center. This type of model has been studied 
in a number of settings, for example analyzing spectral statistics [31 El |35l [36] and the 
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distribution of wavefunctions [SI E]- In general, the results obtained for this particular 
model are more explicit than equivalent results for other graphs so it provides an 
appropriate jumping off point for our investigation. 

A star graph consists of a central vertex of degree B which we refer to as the center 
and B external vertices of degree one which we call the nodes of the star, see Figure [H 
The boundary conditions at the nodes will be either Neumann or Dirichlet respectively, 
where we take the edges to be oriented so x;, = at the node. ?/'^(0) = or ipdif^) = 
where {n}, {d] label the sets of bonds starting at the B^ Neumann or Bd Dirichlet 
nodes respectively. The total number of bonds of the star is then B = B]\f + Bjj and 
we will still use {b} to label the set of all bonds. At the center we choose matching 
conditions that generalize the Neumann boundary conditions: functions are continuous 
at the center, iph{Lb) = for all b, and 

J2^i{L,)=0. (9) 

b 

(j) is not a constant but only convenient notation for the single value of the wave function 
at the central vertex. 




Figure 1. A star (or hydra) graph. 



3.1. The secular equation. 

The Laplace operator on a star graph with Neumann matching conditions at the central 
vertex admits a secular equation, of a particularly simple form, whose solutions form 
the fc-spectrum of the star [20]. From the boundary condition at the nodes and the 
continuity of the wave function at the central vertex, eigenf unctions on the bonds of the 
star have the form 

cos kxn , . sinAjXd 

= '^^^ ' '^'^'''^ = ^m7, ' 

where k'^ is the eigenvalue of the Laplace operator. Substituting in ([9]) produces a secular 
equation, 

tan /cLn — cot = , (11) 
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where the sums are over the bonds starting at Neumann/Dirichlet nodes respectively. 
The positive set of solutions {kj} is the fc-spectrum of the graph. 
If we consider the function 



f{k) = tan kLn — cot kL^ 



f has poles at the points 

■77(771 + 1/2) 



U 



U 77771 

d ^ 



(12) 



(13) 



If the set of bond lengths {Li, . . . ,Lb} is incommensurate, not rationally related, the 
poles are all distinct. As both tan and negative cot are strictly increasing the zeros of /, 
which correspond to eigenvalues of the Laplace operator, are also all distinct each one 
lying between a pair of adjacent poles. Figure [2] shows a schematic representation of /. 




Figure 2. Schematic representation of the functional part of the secular equation 
of the star graph with Neumann matching at the center. 



3.2. The zeta function for Neumann boundary conditions at the nodes. 

As a first step we derive the zeta function of the star where all the nodes have Neumann 
boundary conditions. Let 

/(z) = - VtanzL^ , (14) 
z ^-^ 

n 

where 2 = + it G C. The zeros of / on the positive real axis still correspond to 
square roots of the eigenvalues of the graph Laplacian, while dividing by z removes the 
zero of / at the origin. To formulate our graph zeta function ([1]) we follow a contour 
integral approach introduced in [28], EHl [30] . Each zero kj of / contributes a factor 
to the zeta function and to sum these we use the argument principle |TT] and evaluate 
an integral of the form 

C(^) = ^ [z-'^Q^dz = ^ [z-'^^\ogf{z)dz, (15) 
2vriA f{z) 2771 dz 
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where c encloses the zeros of / and avoids poles which we have already seen are distinct 
when the bond lengths are incommensurate. Figure [3]^a) illustrates the appropriate form 
of the contour c. qj 

(a) (b) 



c 



E](b).Q 



X pole X pole 

o zero d ° zero 



Figure 3. The contours used to evaluate the star graph zeta function, (a) before, and 
(b) after, the contour transformation. 



To analyze C,{s) we deform c to d and integrate along the imaginary axis, see Figure 



Following the contour transformation it is natural to write C,{s) as the sum of two 
terms 

as) = cus) + cp{s) , (16) 

where C/m is the contribution generated by the integral along the imaginary axis and Cp 
is the series over residues arising from the poles of /. At a pole Zq of / we must subtract 
the residue z'^'^'^ . Consequently 

-2s oo 



f $^(m + l/2r 

m=0 



Ch(2.,1/2)5^ f , (17) 



where C,h is the Hurwitz zeta function. 

For the integral along the imaginary axis we get 
1 /"~°° d 

Om{s) = —J (it)-2^-l0g/(it)dt 

t Of course, the contour c has to be thought of as being the limit of a finite contour c„ as n — >■ oo. The 
contour c„ is as in Figure 3(a) but closed by a vertical line with real part a„. The sequence (a„)„gN) 
which goes to infinity, is chosen such that the distance of each a„ to the nearest pole of log/(z) is 
larger than a suitably chosen e > 0. This guarantees that log/(z) remains bounded along the vertical 
lines and so the contributions from the vertical lines vanish for 3f?s>l/2asn— ^oo. The sequence a„ 
can be constructed because of Weyls law, as the average separation of the fc-spectrum is constant. 
I The contour c' is closed by a semicircular arc of radius a„ on the right. Along the arc, f{z) grows at 
most exponentially fast and as n — oo the relevant ratio |/'(z)//(z)| is asymptotically constant and 
contributions from the contour at infinity vanish for 3?s > 1/2. Although we will not stress this again, 
the same construction holds for all zeta functions considered later. 
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fit) 



dt 



where 



f{t) = ^tanhtL, 



(19) 



As t ^ 0, f{t)/t ~ En(^n-^n^V3) + 0(t^) and in this hmit ^ log/(t)/t is proportional 
to t so f lTSj) is vahd at most for 3ft s < 1. Similarly, from the t — )■ oo behavior one obtains 
the restriction 3ft s > 0. 

Our representation of Cim{s) therefore holds in the strip < 3fts < 1. To obtain 
an analytic continuation valid for 3ft s < 1 we split the integral at t = 1 and develop the 
integral over [1, oo), 



c 



Im[S] 



sm TTS 



TT 



dt 



log 



t 



dt+ rr2^^iog(/(t)) dt 



1 

2^ 



(20) 



Collecting these results we obtain the following theorem 



Theorem 4. For the Laplace operator on a star graph with Neumann matching 
conditions at the vertices the zeta function for 3ft s < 1 is given by 



as) = CH{2s,i/2)j2(r 



+ 



n 

smvrs 



-2s 



sm US sm ns 
+ 



27rs 



vr 



-2s d 1 



fit) 



dt 



TT 



oo 1 

-2s ^ 



r^^-iog(/(t)) dt. 



To demonstrate the power of this formulation of the zeta function we calculate 
the derivative of C at zero. We discuss the implication of the result for the spectral 
determinant of quantum graphs in Section [61 Differentiating (120|) we see that 



ClmiO) 



log 



fit) 



t 



+ 



log/(t) 



(21) 



where C/B is the mean bond length. From (fTTl) the pole contribution to C'(0) is given 
by 



CpiO) = 2BCHiO,l/2) + (HiO,l/2)J2 



-2 log 



71 
Ln. 



= -B\og2 , (22) 

where we have used Ci?(0, 1/2) = and 2^^(0, 1/2) = — log(2). Combining the results 
for ^|m(0) and Cp(0), we find 



C'(0) = -log 



2^C 
B 



(23) 
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3.3. Mixed Dirichlet and Neumann conditions at the nodes. 

If we include Bd nodes with Dirichlet boundary conditions the zeta function construction 
can be modified to generate an integral representation of C(s)- Let us define the functions 



/(^) = ^ I XI tan zLn - X cot zL^ j , 



(24) 



f{t) = ^tanhtL^ + ^cothtLd , (25) 

n d 

SO that /(it) = —tf{t). f is defined so that zeros of the secular equation are zeros of / 
but / is not divergent at zero, /(O) = — Yld^~d^ ■ represent the zeta function again 
using the same contour integral 

as) = J^z-'^^dz ^ J^z-'^-^Jogf{z)dz. (26) 

Proceeding as described before, the integral along the imaginary axis is given by 

Cimis) -'^f ^-'^^ log (tm) , (27) 



, . SmTTS 
Clm[S) = 



(28) 



for < < 1. Splitting the integral at t = 1, we obtain the analytical continuation 

-[I i ('^ w) + r ^'w + h 

valid for 3ft s < 1. 

The positive poles of cot zL^ are {mn/ Ld}meN- Summing their contributions as 
well as those corresponding to poles of tan zLn discussed previously 



-2s 



cH.)-E(|:) E('-V2)-+E(^) E.- 

n " m=0 d ^ j=l 

= (2^^ - 1)^(2.) (f ) + Cr{2s) J2 (f ) ■ (29) 
Collecting these results we have the following theorem. 

Theorem 5. For the Laplace operator on a star graph with Neumann matching 
conditions at the central vertex, Bd external nodes with Dirichlet boundary conditions 
and Bn external nodes with Neumann boundary conditions the zeta function for 3fts < 1 
is given by 

C{s) = (2^^ - l)U2s) [jj + Cii(2s) [jj + 



a 



sin7r5 9« d , f ^\ , sinirs 9« d , x , 



Although the addition of bonds with Dirichlet boundary conditions appears a small 
variation of the star graph model - certainly the formulation of the zeta function in 
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theorems H] and are very similar - their addition has a substantial impact on spectral 
properties of the star. This can be seen if we evaluate C'(0)- First we note 



CL(o) 



dt 



log dt + 



1 dt 



log/(t) dt 



log(5^L,i)+log5. 



Differentiating (^p, 



CMO) = 21og2C/?(0)5^ + 2(^^(0)5,, - 2(^(0) log 



TT 



= - log 
Combining the results 

C'(0) = -log 



(30) 



(31) 



^(n^^)(E^.^) • (32) 

In comparison with (l23l) the spectral determinant of a star with Dirichlet and Neumann 
nodes only depends on the lengths of the bonds starting at nodes with Dirichlet boundary 
conditions. 



3.4- Zeta functions with equal bond lengths. 

It will be instructive, for comparison, to evaluate the zeta functions of the star when 
all the bonds have an equal length. In this case the spectrum is known and we can 
calculate the zeta function directly. The results in this case will be more explicit and 
will provide a useful test case when we apply the zeta function to derive the vacuum 
energy of the graph. 

If we consider the star graph where all the nodes have Neumann boundary 
conditions and when the bond lengths are equal, L^, = L for all b, the secular equation 
f llip reduces to B tan kL = 0. The /c-spectrum now consists both of the zeros of 
tan kL, which are simple eigenvalues, and the poles of tan kL, which are eigenvalues 
with multiplicity B — 1. To see that the poles are also eigenvalues one may start 
with the secular equation tan kLb = where the bond lengths are incommensurate. 
Initially every eigenvalue lies between a pair of adjacent poles. As the lengths of the 
bonds are equalized groups of B poles come together trapping B — 1 zeros of the secular 
equation at a pole of tan kL. From this spectrum the zeta function excluding the zero 
mode is 

ji=l m=0 ^ ^ 

= (I) {{B - 1)2^^ -B + 2) Cr{2s) . (33) 
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Consequently 

C'(0) = - 2 log (^) Ci?(0) + (5 - 1)2Ch(0) log 2 + 2C'r{Q) 
= log(^) -(5-l)log2-log27r 

= - log (2^L) . (34) 

This agrees with our previous result with incommensurate bond lengths (!23|) as the 
total graph length is £ = BL. This is expected as the previous result is continuous with 
respect to the bond lengths. 

If we include Bd nodes with Dirichlet boundary conditions keeping all the bond 
lengths equal we can still evaluate the spectrum directly. The secular equation (fTT!) 
reduces to, 

Bn sin^ kL - Bd cos^ kL = . (35) 

If we define a = - arcsin ^Bd/B then zeros of the equation are values of k in the 
set {(mvr ± aTi)/L,m G Z} each element of which is a simple eigenvalue. As in the 
previous case making the bond lengths equal traps eigenvalues at the poles of tan/cL 
and cot kL. The sets {(2m + l)7r/2L, m G Z} and {rmi/L,m G Z} therefore correspond 
to eigenvalues of multiplicity Bn — 1 and Bd — 1 respectively. 

The zeta function can now be calculated directly from the spectrum. We have 

mi=0 m,2=l 

+ (B. - 1) ± a^^^y + (Bo - 1) i: i^y- 

7713=0 7714 = 1 

= (^)~'' [CHi2s,a) + Ch{2s, l-a)+ (^Bd- Bn + {Bm - 1)2'')Cr{2s)\ . (36) 
From this representation of the zeta function it is straightforward to derive C'(0), 
C'(0) = {Bd - 1) log (£j - {Bn - 1) log 2 - {Bd - 1) log(27r) 

+ 21og (^csc(a7r)) - 21og2 , (37) 

where we have used 

CH(0,«) + C/f(0,l-«) = (38) 
CH(0,a) + C^(0,l-«)= log(csc(a7r))-log2. (39) 



From the definition of a we see that csc(a7r) = y B/Bd, consequently 
C'(0) = - 5 log 2 - (fiz) - 1) log(L) + log 5 - log Bd 

._,o,f!™^). (40) 



c 

where C = BL is the total length of the graph. This agrees with (|32|) for incommensurate 
bond lengths if we set Lb = L for all bonds b. 
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4. The zeta function for a star graph with general matching conditions at 
the center 

We introduce the techniques used to study the zeta function of any quantum graphs by 
generalizing the star graph example to a star with Dirichlet conditions at the nodes but 
where the matching conditions at the center have any general form compatible with the 
domain of a self-adjoint Laplace operator. The matching conditions at the central vertex 
will be specified by a pair of B x B matrices A and B using the scheme of Kostrykin 
and Schrader [M], see Section [2l The Laplace operator on the star is then self-adjoint 
if and only if (A, B) has maximal rank and AM^ = BA"!". 



4.1. The secular equation. 

Dirichlet boundary conditions at the nodes imply that an eigenfunction has the form 

ipb{xb) = Cb sin kxb (41) 
on each bond b. The matching condition at the center is defined by the matrix equation 
AV^ + BV' = 0, where i/^ = (^i(Li), . . . , MLb)^ and V' = (-^'i(^i), • • • , -^^(^i?))''- 
Let c = (ci, . . . , cb)^ and define two diagonal B x B matrices 

sin(fcL) = diagjsin kLi, . . . , sin kLs} , (42) 

cos(fcL) = diagjcos kLi, . . . , cos kLs} ■ (43) 
The matching condition can equivalently be written as 

(^A sin(fcL) - kM cos(A;L) j c = . (44) 
k is therefore an eigenvalue if and only if it is a solution of the secular equation 

det (^Asin(fcL) - A;Bcos(A;L)j = . (45) 

We can put this in an alternate form that mirrors the form used for the zeta function 
of a general Laplace operator on a cone [3ll [32] . 

Lemma 1. Let X = diag{a;i, . . . ,xb} and Y = diag{|/i, . . . , i/b} with both X and Y 
invertible. Then 

A 



det ^ ^ j = det (aX-1 - MY'^) \[x,y, 



This is a straightforward consequence of 

x-i 

Y-^ 





(46) 



Applying Lemma [H to the secular equation (H5l) when neither sin kLb or cos kLb is 
zero for any b we obtain 

(a B 

^""^ I (sin(A;L))-i \ (cos(A;L))-i ' " ° ' ^^'^^ 
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Equivalently we may write the secular equation in a form that bears a functional 
similarity to the cases we have already analyzed when the center had Neumann like 
matching conditions, namely 

detf 1 1=0. (48) 

' Ib \ tan(A;L) ' ^ ^ 

4-2. Zeta function calculation. 

Following the final formulation of the secular equation let us define functions 

luu(.L)j- ltanh(tL)J • (^^' 

SO f{t) = /(it). First we note that 

/(O) = /(O) = det ^ j . (50) 

/(O) is generically non-zero in the sense that if /(O) = perturbing the set of bond 
lengths in L by an arbitrarily small amount will make it non-zero. Secondly we will also 
be concerned with the behavior of / in the limit t to infinity. In this limit / has an 
asymptotic expansion of the form 

/(t)~detB + ^ + | + ... . (51) 

We denote by a at the first non-zero coefficient in the expansion, so = if detB ^ 0. 

We are now ready to state the following theorem for the zeta function of the star 
graph with a general matching condition at the central vertex. 

Theorem 6. For the Laplace operator on a star graph with general matching conditions 
at the center, defined by matrices A and B with AB"!" = BA"!", rank(A, B) = B, and 
Dirichlet boundary conditions at the nodes the zeta function on the strip — l/2<3f?s<l 
is given by 

N > / In / 7r \ A^sinvrs sinvrs 9« d , ^ , 

H-^n-4iog(*-/w)d*. 

TT at 



fit) = det 



A B 
Ib j tanh(tL) 

A direct consequence of Theorem [6] is the following simple formula for C'(0), 



B 




C'(0) = -logl— detl ; 1 1. (52) 
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Proof of Theorem Following the zeta function calculations for star graphs with a 
Neumann like matching condition at the center the theorem can be established in a few 
lines. Again we consider the contour integral 

d 



.-4^d. 



-2s 



log f{z)([z 



(53) 



with c as shown in Figure [3] and / defined in (H9l) . We again split the zeta function where 
dm is the contribution of the integral on the imaginary axis and (p is the contribution 
of the poles of / on the positive real axis. 

The poles of / are at (m + l/2)7rL^^ for m integer, the same set of poles as the 
star with Neumann matching at the center and Neumann boundary conditions at the 
nodes, so 



Cp{s) = U2sA/2)Y,(J-^ 



(54) 



On the imaginary axis / (it) = fit) is an even function of t and we have the representation 



Im{S 



sm ns 



TT 



t~'' -log fit) dt. 



(55) 



As t tends to zero fit) ~ /(O) +Cit'^ + . . .. /(O) is generically non zero and consequently, 
for > ( |55l) is valid in the strip < < 1, whereas for A^ = it is valid in the 
strip —1/2 < 3fts < 1. Splitting the integral at t = 1 and subtracting the asymptotic 
behavior of / as t tends to infinity, see flSTl) . we obtain the theorem. □ 



4-3. Comparison with results for a star with Neumann matching at the center. 

For completeness we calculate C'(0) using (152!) when the central vertex has a Neumann 
matching condition in order to establish agreement with our previous results. For 
Neumann like matching condition the matrices A and B can be chosen to be 



/ 1 -1 
1 



A 





-1 



\ 












1 





/ 









1 / 





-1 

Then equation implies ipbiLt) = ipb+iiLb+i) for 6 = 1, ... , B—1 and J2b=i i^bi^b) 
To evaluate C'(0) using fl52|) we need to evaluate /(O) and a^- We first note that 



(56) 



/(O) = det 



det(AL 



(57) 
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Given our choice of A and 



/ Li -L2 
L2 -L3 



/(O) = det 



\ 



: ■•• ■•• ■•• 

... Lb-i -Lb 
V-1 ... -1 -1 -1 / 

/ L2 -L3 ... \ 

L3 -L4 ■•■ : 



(58) 



Li det 





V -1 





Lb-1 —Lb 
-1 -1 -1 / 



B 



Li\{L, 



(59) 



6=1 



Iterating this procedure 

det(AL-B) = LiLadet 

Consequently, we find 

/(O) = det 



/ L3 -L4 
L4 -L5 



\ 





v-1 





Lb~i —Lb 
-1 -1 -1 / 



B 



-L-^^\{L,-L-^^\{L,.m 



6=1 



b=l 



B 



B 



Ib L 



,b=i / \b=i 



(61) 



To obtain the t to infinity behavior of / we rewrite the determinant so 



(62) 



as the rank of A is 5 — 1. As t tends to infinity tanh(i(;L) approaches I^. The first 
non-zero coefficient in the asymptotic expansion of fit) is 



aB-i = det(A-B) = -B 



(63) 



where the determinant can be evaluated following the same expansion used to evaluate 
/(O). Substituting the values of Qb-i and /(O) in f l52|) we again obtain fl32|) . 



5. The zeta function of a general quantum graph. 

A secular equation based on a bond scattering matrix of the graph is a widely employed 
starting point for the derivation of the trace formula of a quantum graph, see e.g. [HI ES] • 
However, the approach we have adopted so far has a natural analogy with techniques 
used to study zeta functions associated to manifolds which we would like to maintain. 
To achieve this it is necessary to employ an equivalent formulation for a secular equation 
associated to a general quantum graph in which the matching conditions appear directly. 
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5.1. Secular equations for general quantum graphs. 

For comparison, the secular equation for a general quantum graph defined in terms of a 
bond scattering matrix (or quantum evolution operator) was introduced by Kottos and 
Smilansky [35| [36] . In our notation the quantization condition on the graph takes the 
form 

det I I - m I L 11=0- (64) 







The scattering matrix S is a unitary 2B x 2B matrix which can be defined in terms of 
the matching conditions on the graph [3^ . 

^{k) = -(A + iA;B)"^(A - iA;B) . (65) 

The following alternative formulations of the secular equation are all equivalent and 
equivalent to (^^; they differ primarily in the form of the functions used to incorporate 
the /c-dependence. We first simply list the secular equations before indicating how to 
derive them. 

The secular equation written in terms of the boundary conditions using exponential 
functions has the form 

/ / -ifcL/2 ifcL/2 \ / -ifcL/2 _ ifcL/2 \ \ 

(A ^i,.,^^ j + ,M j j = , (66) 

Alternatively a similar quantization condition can be obtained with trigonometric 
functions, namely 

det^A^-^^^^^^/^) ^°^(^^/2)V cos(/.L/2) sin(A:L/2) \\ _ 

sin(A;L/2) cos(A;L/2) J + - cos(A;L/2) sin(A;L/2) JJ" ^ ' 

The final formulation that we employ later has a similarity with the secular equation of 
the Neumann star graph (fTTl) . 



Derivation of (!67|) : the matching conditions on the graph are specified by the matrix 
equation 

Ail) + Mip' = . (69) 

On each bond the eigenfunction has the form 

i^bixb) = Cb sin kxb + Cb cos kxb ■ (70) 

Consequently the vectors of values of the function and its derivative at the ends of the 
intervals corresponding to each bond are 

= {ci, ... ,cb,ci sin kLi + ci cos kLi, ... ,cb sin kLs + Cb cos kLs)'^ , (71) 

jp' = k{ci, . . . , Cb, ci sin kLi — c\ cos kL\, . . . ,CBsinkLB — cb cos kLB)^ ■ (72) 
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Writing c = (ci, . . . , cb)^ and c = (ci, . . . cbY the matching conditions (1^ take the 
form 

^ ( sin(A;L) cos(fcL) ) ( c ) + ^® ( - cos(A;L) sin(A;L) ) ( c ) = ° ' ^^^^ 
Eigenfunctions therefore exist for values of k that solve the secular equation, 

^ ^ sin(A;L) cos(A;L) ^ ~'~ ^ — cos(A;L) sin(A;L) ^ 

This can be reduced to the form given in (!67|) post-multiplying by 

cos(A;L/2) sin(A;L/2) 
-sin(A;L/2) cos(A;L/2) 

a matrix with unit determinant. The secular equation ( l66l) is obtained from ( l67l) after 
post-multiplying by 




Alternatively the same result comes from expressing the wavefunctions on the bonds as 
sums of plane waves. Finally equation ( 168|) can be obtained from (17^ by multiplication 
on the right by 



Is ^ ^_/-cot(A;L) csc(A;L) 

sin(A;L) cos(A;L) ) ~ 1 1^ 



(75) 



This matrix has a nonzero determinant and the inverse exists provided k is not a zero 
of sin kLi, for any h. 

5.2. General zeta function calculation. 

To derive the zeta function it is convenient to start with the secular equation written in 
the form ( |68|) which resembles the equation ( ITTj) for a star graph that we began with. 
In the general graph case we define the two functions 

/W^det|A-.B| —7 II. (76) 





such that f{t) = /(it). Poles of / lie at integer multiples of n/Lb on the real axis. When 
the set of bond lengths {Lh} is incommensurate the poles of cot zLh are all distinct. 
However, as there are four functions with the same set of poles in the determinant it 
is possible for a linear combination of the functions to cancel some of the poles of cot. 
Generically, however, the set of poles is {rmr / Li,\m G Z, 6 = 1, . . . , B} and this is the 
case we consider. 
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Again we represent using the contour c, see Figure El 



C(^) = / ^ '7(iy = j ^~''j-^^ogf{z)dz. (78) 



Notice, that at the origin 



/(0) = det|A-B| ||=/(0). (79) 



L 1 -L 1 
L 1 L 

If /(O) = an arbitrarily small change in the bond lengths can make this non-zero, so 
the integral formulation naturally excludes zero modes. 

As before, we split the zeta function into the sum of the imaginary axis integral 
and the contribution of the poles on the real axis C,{s) = Cim{s) + Cp{s). As we have 
previously seen in the case of the star graph for this set of poles 



B 

\ ( TT 



6=1 



-2s 



Cp(5) = Ci?(2s) > ( - . (80) 



On the imaginary axis 



sinvrs /"°° ._2s d 



C.Ub) = / t-^^-log/(t)dt . (81) 

TT Jo 

As t tends to infinity, up to exponentially damped terms, 

f(t) ~ det (A - tl) = det B t^^ + c^B-it^^'^ + . . . + Cit + det A (82) 

Let ctv be the first non-zero coefficient of the highest power of t, i.e. = 2B when 
det B 7^ 0. Splitting the integral at t = 1 and subtracting the t to infinity behavior we 
obtain the following integral formulation of the zeta function. 

Theorem 7. For the Laplace operator on a graph whose vertex matching conditions are 
defined by a pair of matrices A and M, with AB"'' = BA"'' and rank(A, B) = B, the zeta 
function on the strip — 1/2 < 3fJs < 1 is given by 

N > / N / TT \-2« A^sinvrs sinvrs 9c d , ?/ x , 

c(.) = c«(2.) E (^) + + — f t--- log/w it 



fit) = det 



coth(tL) — csch(tL) 
— csch(tL) coth(tL) 

where the poles of f{z) = f{—iz) are the whole of the set {rrnr / Lf,\m G Z, 6 = 1, . . . , B}. 

Theorem [7] implies the following simple formula for the derivative of zeta at zero in 
terms of the matching conditions at the graph vertices. 



C{0) = -\og 



CAf nf=l 



det A 




(83) 



We have now obtained results for the spectral zeta function of a generic quantum 
graph along with a number of specific zeta function calculations for the formative case 
of a star. In the balance we assess the consequences of these results for the spectral 
determinant, vacuum energy and heat kernel asymptotics of quantum graphs. 
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6. Comparison with results for the spectral determinant. 



The spectral determinant of a Schrodinger operator on a quantum graph is an important 
subject in its own right which has previously been investigated by a number of authors 
, [17] . Formally the spectral determinant of the graph Laplacian is 



det'(-Al 



j=0 



3 ■ 



The zeta function representation presented here allows a direct evaluation of the 
regularized spectral determinant, det'(— A) = exp(— (^'(0)). Consequently, ( l83l) provides 
the spectral determinant of the Laplace operator on a general graph in terms of the 
vertex matching conditions presented in Theorem [H 

We wish to compare our results derived from the zeta function with the literature. 
As the known formulations of the spectral determinant each apply to certain classes 
of quantum graphs for the purpose of comparison we take the example of the Laplace 
operator on a star with Neumann conditions at all the vertices. From the star graph 
zeta function, see (l23il . in this case 



det'(-A) 



B 



5) 



Friedlanders formulation for the spectral determinant of Schrodinger operators on 
graphs with delta type vertex matching conditions [T7] reduces, in the case of Neummann 
like conditions at the vertices, to the following spectral determinant. 



detV-Al 



11, 



det'R. 



6) 



The notation det' denotes the determinant excluding eigenvalues of zero and d{v) is the 
degree of the vertex v. R is the Dirichlet-to-Neumann operator at zero energy, slVxV 
matrix defined by 

V w 

V = w 



(87) 



When V ^ w the sum in (IHTj) is over bonds connecting v and w and when v = w one 
sums over bonds originating at v excluding loops. For the star graph V = B + 1 and 
the matrix R can be written in the form 



R 








\ 



L 



B 
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It is convenient to combine Y\.b with det'i? so 



B 



Yl Lb det'R = det' 



v 



B j-1 




... 




-1 


1 








-1 





1 





-1 








1 



(89) 



Applying column operations it is clear that this matrix has eigenvalues and 1 with 
multiplicity 1 and B respectively. Consequently, for the star graph (fl^ Lb) det' R = 1 
and 

det'(-A) = . (90) 



This does not agree with flSSl) . the difference between the results being the factor V. 
Restricting to the interval with Neumann boundary conditions, the associated zeta 
function is 

(is) = {iy'\R{2s) (91) 



and 



(92) 



C'(0) = -log(2L). 

This situation corresponds to B = 1, C = L and V = 2. Consequently (!85|) is seen to be 
correct differing from (!90l) by some normalization constant. Finally, it is worth noting 
that Theorem [1] agrees with the conjecture of Texier for the spectral determinant of a 
quantum graph with general boundary conditions 



7. Vacuum energy and Casimir force of graphs. 

The vacuum (Casimir) energy of quantum graphs has also been a topic of recent research 
[H |19] . A quantum graph provides a simple model of vacuum energy for a system whose 
corresponding classical dynamics is ergodic. The quantum graph model also naturally 
introduces a set of independent length scales, the set of bond lengths, so techniques 
applicable to quantum graphs may provide insight into the approaches required to 
tackle less symmetric but physically relevant structures in thee dimensions. However, the 
quantum graph vacuum energy should only be thought of as a mathematical construction 
rather than a lower dimensional model for the vacuum energy fluctuations of some 
underlying physical configuration of narrow tubes. It remains to be established how the 
vacuum energy of a graph is related to the vacuum energy of a thickened manifold with 
specific boundary conditions built over the graph skeleton. 

Formally, the vacuum energy associated with the Laplace operator on a graph is 
given by 

^ oo 
j=0 
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The zeta function regularization of this sum is = ^({—1/2). However, not all the 
zeta function expressions constructed so far are valid for 3ft s = —1/2 and these cases 
will need to be examined with more care. As we will see, C(— 1/2) generically diverges, 
but the Casimir force acting on a bond will be well defined. 

The vacuum energy of a star graph with equal bond lengths was investigated by 
Fulling, Kaplan and Wilson in [19] using the heat kernel regularization. Of particular 
interest they discovered that the sign of the Casimir force depends on the number of 
bonds and for star graphs with more than three bonds the force is repulsive. In |1] 
Berkolaiko, Harrison and Wilson applied the quantum graph trace formula to obtain 
periodic orbit expansions of the vacuum energy of general graphs. The zeta function 
approach is new and complimentary, relating the graph vacuum energy directly to the 
matching conditions at the graph vertices. 

7.1. Star graphs with equal bond lengths. 

We begin with the simple case of a star with Neumann matching conditions at the 
center and Neumann boundary conditions at the nodes. The zeta function for equal 
bond lengths L was given in equation fl33l) . Evaluating the vacuum energy 

= ^{3 - B)U-l) 

^ (5-3) (94) 



48L 

which agrees with the results in [U |T9]. The Casimir force is proportional to the 
derivative of Ec with respect to the bond length L and we see that for 5 > 3 
increasing the bond lengths is energetically favorable and consequently the Casimir 
force is repulsive. 

If we now allow Dirichlet boundary conditions at some nodes, ({s) was given in 
equation fl36l) . and we find 



2(C//(-l, a) + Ch{-1, 1 - «)) + {2Bd -Bm- 1)C,r{-1] 
■{B^ - 2Bn + 1) + ^(CHi-l, «) + Cni-l, 1 - «))• (95) 



E.= - 
' 4L 

TT 

Using the relationship between the Hurwitz zeta function at negative integers and 
Bernoulli polynomials one can obtain, 

Ch{-1, a) + Ch{-1, 1 - a) = + a - a^ (96) 



from which we can write the vacuum energy in a concise form as a function of B and 
Bo, the number of nodes with Dirichlet boundary conditions, 

^c=^{b- HBn + 1) + 24a(l - a)) . (97) 

Equation f l97|l generalizes the previous results for star graphs. For example on the 
graph with three bonds the vacuum energy is positive if -Bd = 1 and negative for 
Bo = 2. In general increasing the number of Dirichlet nodes changes the vacuum energy 
from positive to negative and consequently changes the Casimir force from repulsive to 
attractive. 
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7.2. Star graphs with incommensurate hand lengths. 

For star graphs where the set of bond lengths is incommensurate the presentations of 
the zeta function provide an integral formulation of the vacuum energy in contrast to the 
conditionally convergent periodic orbit sum obtained by Berkolaiko et al. [4J . For a star 
with a Neumann matching condition at the center and mixed Dirichlet and Neumann 
boundary conditions at the nodes from Theorem [5] we obtain the vacuum energy 



t 



Y,n tanh tLn + Y.d ^oth tLa 



■dt .(98) 



A similar formula holds when all the nodes have Neumann boundary conditions. 
For comparison, if we set the bond lengths equal, (p8|) reduces to 
TT /„ \ L BNsedi^itL) - BDCsch\tL) 



E. 



B-Wd 



dt . 



(99) 



48LV " 27: J Q " Eat tanh(tL) + 5^ coth(tL) 

This is an integral formulation of the vacuum energy evaluated previously (1971) . As an 
indirect consequence we have obtained the following integral 

2 



a sech x — h csch x 

X ; — ax 

a tanh x + b coth x 



sm 



sm 



-1 



a + h 



TT 



(100) 



This does not appear to have been known. 



7.3. Casimir force on a star with a general matching condition at the center. 

The integral formulation of the graph zeta function obtained in Section HI was for values 
of s in the strip — l/2<3fJs<l. To analyze the vacuum energy we must first continue 
the zeta function formula to include s = —1/2. The restriction —1/2 < 3fJs came from 
the behavior of the function / at infinity. 

A 

t^ 



fit) = det 



Ib \ tanh(tL) 



(101) 



where was the first nonzero term in the asymptotic expansion of / and a^^j is the 
second nonzero coefficient, generically j = 1. From this we find 



log/(t) ~ log 



+ ^ + o(r(^-+i)) 



Recall 



C/m(s) 



smvrs 



-2s 



dt 



log fit) dt 



(102) 



(103) 



When we develop the imaginary axis integral subtracting the leading and subleading 
order behavior at infinity leaves an integral convergent for 3? s > — (j ' + 1)/2 as required. 
We arrive at the zeta function formula 



B 



as) = CH[2s 



TT 



6=1 



2s 



+ 



sm TTS 



TT 



t' 



-2s 



dt 



log fit) dt 



+ 1 't''4(^og(t^fit))-^)dt-^- , , . 

/i dt V ^ ^ «JvtV 2s a7v(2s+j)_ 



(104) 
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The vacuum energy E^^ is divergent if j = 1 which would be the general case. However, 
as terms dependent of the bond lengths are exponentially suppressed in the t to infinity 
behavior both and aN+j are independent of the set of bond lengths {Lh\. Changes in 
the vacuum energy and the corresponding Casimir forces are the observable quantities. 
To find the Casimir force on bond (3 the vacuum energy is differentiated with respect to 
the bond length Lp and therefore the Casimir force is well defined, namely 

It is simple to differentiate / with respect to as only one term in the matrix 
determinant depends on the length of /3. 

Although the form of the Casimir force looks somewhat different from our previous 
results (l97|l and (198|) we can compare them in the case when all the bond lengths are 
equal. Inserting the Neumann matching matrices (13^ . we observe 

1 ^ 

/W = -^E n tanhtL,-. (106) 

6=1 j=B,j^b 

Consequently N = B — 1 and a^+j = 0. Differentiating with respect to and setting 
Lb = L for all b E B we find 

9 = (107) 



dL0 ' B tanh(tL) 

Substituting in fllOSp and integrating we obtain a Casimir force 



TT 



,-(3-25). (108) 

The vacuum energy of a star with a Neumann center and Dirichlet nodes f l97|) when 
Bd = B and consequently a = vr"^ arcsin ^Bd/B = 1/2 is = 7r(3 - 2B)/A8L. 
Differentiating, the Casimir force agrees with fllOSp which is the force on the bond P 
rather than the force on all B bonds of the star. 



7.4- Casimir force on a general quantum graph. 

We conclude with the equivalent result for the Casimir force on a bond of a general 
graph in terms of the matching conditions defined on the graph. Here again we first 
continue the zeta function formula to a form valid at s = —1/2. The vacuum energy will 
still generically be divergent as is the case with general matching conditions in a star. 
However, the asymptotic behavior of / is independent of the bond lengths as there is 
no potential or curvature on the graph and consequently divergent terms do not appear 
in the Casimir force. 

To obtain a form of the zeta function valid for s = —1/2 we again subtract both 
the leading and subleading order behavior of f{t) in the integral along the imaginary 
axis, which follow from 

coth(tL) — csch(tL) 
— csch(tL) coth(tL) 



/(t) = det A - tB 
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~ det(A - tn) = CNt^ + cn-j r-^ + o(r-^-^), (109) 

where cat and c^^j are the highest order non-zero coefficients in the expansion of 
det(A — tM). From here, we continue noting 



iog/(t) ~ log(c^t^) + ^ + 0(r(^+^)) 



CnP 



(110) 



Subtracting the leading and subleading order asymptotic behavior from / in (11031) leaves 
an integral convergent for 3f?s>— (j + l)/2. The zeta function has the form 

B 



(is) 



2s 



+ 1 t-'-i-A^ogit-'^fit))- 



smvrs 



TT 



t' 



-2s 



-N 



dt 



log fit) dt 



CN-j J 



dt-\ , , . 

2s cn{2s + j) 



(111) 



dt 

In general j = 1 and Ec = ^C{—l/2) is divergent. The coefficients c„ are independent 
of the bond lengths. Differentiating with respect to Lp and setting s = — 1/2 we obtain 
the Casimir force on the bond /3 stated in Theorem [2l 



8. Heat kernel asymptotics. 

The heat kernel K{t) has already been mentioned in the context of an alternative 
renormalization scheme for the vacuum energy. It is defined by 

oo 

K{t) = Y,e~^'' (112) 
i=i 

and in the one-dimensional setting considered here it is known to have an asymptotic 
expansion as t — )■ of the form 

oo 

K{t)^ eet'-'/\ (113) 

£=0,1/2,1,... 

The heat kernel coefficients ae are related to the zeta function through the following 
connection ^2j, 

ee = ResiC{s)T{s))U=i/2-i ■ (114) 

If we consider the heat kernel of a star graph with a Neumann matching condition 
at the center the zeta function representation can be used to evaluate the heat kernel 
coefficients. We find 



c 1 



Neuman nodes, 



K(t) ~ < ^^'^^ (115) 
I ^= — ^(Bf) — 1) Dirichlet and Neuman nodes. 



2 

C 

These results are independent of the incommensurability of the bond lengths. 

When considering the asymptotic t — )■ expansion of the heat kernel for the other 
cases we have to realize that the representations of the zeta functions are only valid for 
—1/2 < < 1; see Theorems and [71 In order to find the heat kernel coefficients for 
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the general case we have to provide an analytic continuation that is valid further to the 
left of that strip. 

In detail, for the star graph with general matching conditions at the center, as 
shown previously, we have 



f{t) = det 



with 

A B 
Ib \ tanh(tL) 

The restriction —1/2 < comes from the last integral and it is the t — t- oo behavior 
of /(t), namely 

/(t)~detB + ^ + | + ..., (116) 

that produces the restriction. By adding and subtracting the t — )■ oo behavior of the 
integrand in the last integral the analytical continuation can be found. 
Consider the case detB ^ first, then we have the structure 

log/(t)~log(detB) + 5^-J, (117) 

n=l 

which defines the numerical multipliers 6„ using equation flll6p . So furthermore 

i^o,m~-±%- (118) 

n=l 

Adding and subtracting the leading asymptotic terms of this expansion we write 

C,).C«(24)i:(^)-%^/r^-l.o,/Wd. (n. 



oo 



Sm TTS 



TT 







n=l 



sinfvrs) nhn 



dt ^ — ^ > 



TT ^ 2s + n 



n=l 



a representation valid for — (A^ + l)/2 < JJs < 1. This is the form we use, together 
with equation (I114p . to find the heat kernel coefficients. For s = 1/2 only the first term 
contributes to the residue and thus oq. All other heat kernel coefficients are determined, 
by construction, by the last term. 
From equation flllQp we read off 

Res C (i) = ^. C(0) = 0. 



C(-n) = (-l)"+in62n, neN. 
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For the heat kernel coefficients this shows 

2a/ TT 



52^+1 02£+l 

TT 



&2n 



■5 r(n) 

The full asymptotic heat kernel expansion therefore reads 

Once the matching conditions are fixed the numerical coefficients 62fc-i are easily found 
from equations f lll6p and f lll7p using an algebraic computer program. 

If detB = 0, assume is the ffist non- vanishing coefficient in equation flll6p . The 
relevant expansions in that case are 

00 

fit) ~ 

n=N 

(00 
1 + V ^t- 
00 

= - A^logt + loga^ + ^M"", 

n=l 

where the numerical multipliers bn are defined by the last equation. Therefore, 

n=l 

Comparing with equation f lll7p the only difference is the additional first term which 
changes ({0) = ei/2 = into ({0) = —N/2. The heat kernel expansion now reads 

so that = reduces to equation fll20p as it must. As before, once the matching 
conditions are fixed the numbers 62fc-i are easily found. 

Exactly the same calculation goes through for the general graph with general 
boundary conditions. The relevant t — )■ 00 behavior this time is given by equation 
( l82l) . Let denote the highest power in this expansion; the structure then is 

N N 

CN-e 



m -yc^e = t^Y^ 



n=0 e=o 

00 



log fit) N\ogt + \ogCN + y^., 



n=l 
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where the last equation defines the multiphers 6„. This shows 

^log/(t)~--)^^, 

71=1 

and from the previous calculation we can immediately write down the heat kernel 
expansion 

9. The piston graph 

To demonstrate how the techniques introduced in this article come together in a specific 
example of interest we apply the previous results to the case of a star graph with two 
bonds where there is a delta-type matching at the central vertex and Dirichlet conditions 
at the nodes. The first bond has length Li = L and the remaining bond has length 
L2 = C — L. One can therefore consider the graph as related to a narrow piston in the 
limit that the width shrinks to zero, see Figure HI Certainly in this limit the spectrum 
of an epsilon thick piston with appropriate matching condition on the piston membrane 
can be made to approach that of the corresponding graph. In general quantum graph 
vertex matching conditions are approximable through a suitable choice of a Schrodinger 
operator on an epsilon thick neighborhood [16j. Consequently we refer to this example 
as a piston graph where we regard the vacuum energy, and related Casimir force, as 
a function of the bond length L keeping the total length of the piston C fixed. The 
piston membrane, the central vertex, may then move on the graph approaching one of 
the nodes keeping C constant. The piston model of vacuum energy was introduced by 
Cavalcanti [10], for recent developments see [T8l[22ll23l[271lin Ii5] . 



(a) 



e 

r 



L C-L 

(b) 



Figure 4. (a) An e thick piston, (b) The corresponding piston graph. 



Specifically we consider the negative Laplace operator on the piston graph where 
■?/'i(0) = '^2(0) = and at the central vertex 'ipi^L) = ip2{C — L) = cj) and 
ip'i{L) + ip'2{C — L) = X(f). Again is simply convenient notation to denote the value of 
the continuous graph function at the central vertex. These boundary conditions define a 
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self-adjoint operator on the graph. Following Section H] we write the matching conditions 
at the center using the A, B matrix pair, namely 



1 -1 
-A 




(124) 



Such matching conditions produce an energy dependent unitary scattering matrix for 
the central vertex, 

s{k)= -{K + imy^ik-im) (125) 

1 [ \ 2ik 



2ik-\ \ 2ik A 



(126) 



The transformed secular function f{t) was given in terms of the matching conditions 
conditions (H^ . 

- tanh Lt tanh(£ — L)t 

t t ^ ' t^ 

We see ai = —2 and 02 = —A as f(t) ~ —2t~^ — Xt~^ in the large t limit. Substituting 
into the zeta function representation we get 

'\ d 



-- tanh Lt - - tanh(£ - L)t - — tanh Lt tanh(£ - L)t. (128) 



C(.)=CH(2.,^)7r-^^(L^^ + (£-L)^^)+^ / t-^^^log/(t)dt 

for - 1 < < 1. (129) 



^ log /(t) dt . (130) 



dt \ \ 2 J 2t 2s 2{2s 

Evaluating the Casimir force we will differentiate with respect to L keeping the total 
length of the graph C fixed. The Casimir force is then 
7r£(2L - C) J_ d 
"~ 48L2(£-L)2 ^2^ 

The Casimir force on the central vertex of the Piston is plotted for some representative 
choices of the coupling parameter A in Figure [5l A positive force Fc acts to increase L so 
we see that the central vertex is attracted to the nearest node for all coupling constants A. 
Increasing A increases the magnitude of the attraction without fundamentally changing 
its nature. The effect rapidly saturates for values of A above A = 100. 

Note that the result for the piston graph does not contradict that for the Casimir 
force on a star graph with equal bond lengths. As we have fixed the total length 
of the piston the force when the bond lengths are equal is clearly zero by symmetry. 
Nevertheless the force on a star with two bonds, constrained to have equal lengths, is 
still attractive. 

The asymptotic expansion of the heat kernel of the piston graph is determined by 
the t — )■ 00 behavior of f{t) ~ — 2/t — A/t^. Applying fll22p in this case the asymptotic 
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Figure 5. Casimir force on on the central vertex of the piston graph plotted for some 
representative choices of the coupling parameter A. The total length of the piston 
£ = 1. When A = the derivatives on either side of the central vertex match and the 
vertex is invisible. In the limit A — > cx) the graph decouples into two intervals with 
Dirichlet boundary conditions, = 'kH{2L - C)I2AL'^{C - Lf. 



expansion of the heat kernel as t — )■ reads 

C 1 

jA^t 2 



E 

fc=l,3/2,2,. 



-1 



\2k 



{2k-l)T{k-\) \2 



A 



2k-l 



j.k-1/2 



(131) 



10. Conclusions 



In this article we have developed a systematic contour integral technique to analyze 
spectral zeta functions on graphs. The main strength of the formalism, apart from 
its simplicity, is that one can analyze general graphs with general vertex matching 
conditions in one calculation. Particular graphs and boundary conditions are extracted 
with ease. 

The secular equation, together with the argument principle, are at the center of the 
approach. They allow to write down an integral representation for the zeta function valid 
in the half-plane 3fJs > 1/2, see equations ( 1761) and ( 1781) . The analytic continuation to the 
whole complex plane is obtained from an asymptotic behavior of the secular equation, 
see Section 8. As a result it is straightforward to evaluate the zeta determinant, equation 
(1831) . the Casimir force. Theorem [2], and the heat kernel coefficients, equation (I123p . for 
which we obtain new results in greater generality than those in the literature. 
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